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A FISSION GAS RELEASE CORRELATION FOR URANIUM NITRIDE FUEL PINS 
by Michael B. Weinstein and Harry W. Davison 
Lewis Research Center 

SUMMARY 

A mathematical model has been developed to predict fission gas release rates from 
uranium nitride (UN) fuel pins clad with various materials. Predicted and measured 
gas releases have been correlated for UN pins tested to 4. 6 -atom -percent fissions 
(burnup) at temperatures which primarily ranged from 1250 to 1660 K. Most of the fuel 
pins were clad with Nb-lZr or PWC-11; however, several were clad with W-25Re or 
T-111. 

The model assumes gas release from individual grains to the grain boundaries. 

For grains at the fuel surface, gas release is instantaneous. For grains inside the 
fuel, the gas collects at the grain boundaries where the gas volume grows and eventually 
interconnects with the outside surface of the fuel. The fraction of grain bovmdary poros- 
ity which has interconnected has been assumed to be proportional to the porosity, with 
the proportionality constant calculated from the gas release data. As a similifying as- 
sumption it is assumed that the pressxire of the intergranular gas is equal to the fuel- 
clad contact pressure. Fuel and clad strength properties were used to calculate the 
contact pressvire. Empirically, interconnection only occurs if the fuel is above a set of 
minimum temperature or bxirnup conditions. Diffusion coefficients for the intergranular 
transport of fission gas were calculated from fission gas release rate data measured in 
a sweep gas facility. 


INTRODUCTION 

Reliable methods to predict the amount of fission gas release are required for the 
most efficient design of space power reactor fuel elements. A program to develop these 
methods for xiranium nitride (UN) fuel pins has been in progress at NASA Lewis 
Research Center. 

Model development is difficult due to the complex and interacting fission gas trans- 
port and release processes operating in fissioning fuel elements. Carroll, in a recent 



review (ref. 1), cites a multiplicity of different transport, trapping, and release mech- 
anisms and models which have been used to explain fission gas behavior in fuels. Pos- 
sible fission gas transport mechanisms include atomic diffusion, fission-enhanced dif- 
fusion, and gas bubble diffusion, all of which can be complicated by trappii^ and release 
processes. 

Atoms diffusing through the fuel can be trapped by defects or by stationary gas bub- 
bles. They can be released from these traps by defect annealing or by kinetic resolu- 
tion from the bubbles. Similarly, bubbles diffusing under the influence of a driving 
force, such as a temperature gradient, may be trapped or delayed by solid fission 
product precipitates, dislocations, or grain boundaries. Release from these traps oc- 
curs when the driving force exceeds the retarding force. 

These transport-trapping processes are themselves the product of multiple causes. 
As an example, Carroll (ref. 1) points out that possible bubble motion mechanisms are 
surface diffusion, volume (or bulk) diffusion, and evaporation-condensation. Domina- 
tion of one bubble diffusion mechanism over the others depends on the material proper- 
ties, the fuel temperature, and the bubble size (Nichols, ref. 2). 

Fission gas release mechanisms include direct recoil (the ejection of fission- 
produced atoms from the fuel) , knockout (ref. 3) (the ejection of gas and fuel atoms from 
the fuel surface by recoiling fission fragments) , as well as the release when diffusing 
gas atoms or bubbles reach a free surface. Free surfaces include internal cracks or 
grain boundaries open to the exterior surface of the fuel. 

Tests to study the variation of radioactive fission gas release rates from UN with 
three major test variables - temperature, fission-rate density, and burnup - were con- 
ducted usii^ a sweep gas facility located at the NASA Plum Brook Reactor Facility. A 
description of the sweep gas facility has been given in reference 4. 

In these experiments, small-size, high -enrichment UN samples were irradiated at 
high fission-rate densities at temperatures from 425 to 1780 K. Based on a preliminary 
diffusion transport model (atomic plus fission -enhanced diffusion with desorption and re- 
coil release), equations were developed to relate measured fission gas release rates to 
the test variables. Unknown constants in these equations were evaluated from the data 
(ref. 5). 

With the completion of additional tests and analyses, transport by bubble diffusion 
was added to the preliminary model, and the agreement between measured and calculated 
release rates improved. 

In actual fuel pins, the restrainii^ effect of the fuel element clad has a significant 
effect on fuel swelling (ref, 6) and would also be expected to significantly affect fission 
gas release. CYGRO-3 (ref. 7) a new fuel rod performance program includes this clad 
restraint effect on gas release by linking gas release to bubble size which is in turn de- 
pendent on fuel element stresses. 
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Fission gas release from fuels has been linked to the interconnection of grain bound- 
ary porosity. Carroll (ref. 1) mentions the interconnection concept as is used to explain 
UN behavior, while Ainscough (ref. 8) pictorially shows the development of intergranular 
bubbles and cracks in UO 2 . 

We have developed a model relating gas release to the amount of interconnected 
grain boundary porosity, the growth of which depends on the fuel stress as determined by 
the fuel-clad contact pressure. 


RESULTS AND DISCUSSION 

Equations have been developed and correlated with measured fission gas releases 
from UN fuel pins representing various fuel -clad combinations. The development of 
these equations can be summarized as follows; 

First the release rates of radioactive fission-product noble gas isotopes (krypton 
and xenon) from four UN fuel specimens were measured using a helium sweep gas 
facility, a description of which is given in reference 4. (Test operations, data -taking 
procedures, preliminary analyses and observations are discussed in ref. 5. Pertinent 
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sweep gas data is presented in appendix A.) The measured krypton-88 (Kr ) release 
rates were then correlated with a gas release equation which was derived by considering 
diffusional transport and release from fuel specimens. Three separate diffusion mech- 
anisms are included in the gas release model: fission-enhanced diffusion, atomic dif- 
fusion, and bubble "diffusion” (bubble motion caused by surface diffusion). 

The postirradiation photomicrographs of our irradiated fuel specimens show exten- 
sive interconnection of intergranular porosity for the specimens operated at ~1500 K or 
higher. Hence, we linked gas release from these specimens to the gas release from in- 
dividual grains to the grain boundaries. 

The model which was then developed to represent the release of stable (nonradio- 
active) gas from clad fuel pins considers gas (at a constant pressiire) being released to 
the individual grain boundaries where the gas volume increases and eventually inter- 
connects to the outside. As simplifications we assume all the fuel porosity to be located 
at the grain boundaries, and assume the amovint of interconnected porosity to be propor- 
tional to the porosity. We also assume that the pressure of the intergranular gas is 
equal to the fuel -clad contact pressure which can be approximately calculated from fuel 
and clad creep strength and other material properties. The stable gas release model is 
then correlated with data obtained from various sources. 

Krypton-88 is used for most of oim experimental analysis because its 2. 4-MeV 
gamma ray is easiest to monitor and coimt. 
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The detailed derivations of all the correlation equations are carried out in appen- 
dixes B, C, and D. The stable gas release data from outside sources is listed in appen- 
dix E. 


Radioactive Fission Gas Release Rate Correlation 
88 

In an earlier report (ref. 5) Kr release rates from three UN specimens were cor- 
related with the predictions of a preliminary fission gas transport -release model which 
included the mechanisms of atomic diffusion, fission-enhanced diffusion, and direct re- 
coil release. Postirradiation photomicrographs of the fuel specimens indicated that 
bubble migration was also a probable gas migration mechanism, and that gas transport 
was first to the grain boundaries (fig. 1), then out of the pin. Additional gas release 
data from a fourth UN specimen (capsule 124) reinforced the earlier conclusions. Data 



Outside Inside 

radius radius 

Figure 1. - Section of fuei pin in capsule 123. X500. 
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TABLE I. - OUTLINE OF IRRADIATION TESTS 


Capsule 

Total ir- 
radiation 
time, 
hr 

Total 

burnup, 

percent 

Primary 
test tem- 
perature, 
K 

Temperature range 
(side thermocouple), 
K 

Fission-rate density 
range, 
fissions 
cm^-sec 

Number of 
sweep gas 
samples 

121 

1294 

8. 3 

1780 

1215 to 1789 

3.34 to 8.29xlo‘‘^ 

9 

122 

1923 

6.6 

1223 

593 to 1225 

.71 to 3.64 

19 

123 

1894 

6.0 

1505 

579 to 1512 

.49 to 4.45 

19 

124 

1751 

7. 8 

Several 

^425 to 1323 

.36 to 7.57 

15 





‘’525 to 1357 

.35 to 3.15 

11 


3 , 

Helium sweep gas. 
'^Neon sweep gas. 


for all four capsule tests is presented in table I and in detail in appendix A , 

The derivation of the following fission gas release rate equation which includes a 
term for bubble motion is presented in appendix B: 



= (1 + Qf^b) 


a,YSLf 
^ 1 fYS 

( -Qq/RT -Q_/Rt\ 

Df J + D e ° ^ e ® 

10 0 0 0 

K“T^ 

1/2' 

4T 

X 

\ 



\ / 



where b > 1. 5 percent. (Symbols are defined in appendix F.) 
The bubble motion term in equation (1) which is 


K^T^ 


-Q^/RT 
e ® 


is an approximation which neglects bubble growth, size dependent velocity, and gas 
transport to and from the matrix. The term does include power, thermal conductivity, 
and temperature dependences. 

2 2 2 

Equation (1) contains seven unknown constants: a^, a^, ’ ^o’ ^2® ’ 

and Q . The values for these constants were obtained ty conducting a least -squares 
s 



fitting analysis on equation (1) and the Kr®®-helium release rate data given in appen- 
dix A. Steps outlining the use of such a least-squares analysis with a nonlinear equa- 
tion,, such as equation (1), are presented in appendix C. 

The resulting gas release equation is 


Rpin = (1 + 0-376b) 


18 800 


3. 6 5x10 "^Yf 


fY 

.1/2 


1.55xl0“^®f + 4. 49x10 "®e 


+ 


1.89xl0~^'^f 

k2t2 



(2) 


where b > 1.5 percent, ^2 “ ® capsules 122 and 124. Equation (2) is given next in 

OO 

terms of the activity of Kr°° in the sweep gas. 


Aj^ = (1 + 0.376b) 


3.0 7X10 ~^f 

4T 


18 800 


+ f |1.59xl0“^^f + 4.61X10 


.-15, 


+ 


1.94xl0'23f 

r2t2 



(3) 


where = Rp^^(Kr® V4. 285x10^, x = 6. 877x10"^ second'^, and Y = 3.6xl0"2. The 

Q QQ 

constant 4.285x10 arises from the decay of Kr in traveling from the fuel specimen to 
the detector. The calculated activities are compared graphically to our data in figure 2. 
For convenience, calculated values for all of the data points are included in the appro- 
priate tables of appendix A . 


6 




For the data points included in figtire 2, the mean absolute relative difference, 





. is 13. 8 percent, and the standard deviation is 11.1 percent. 

As noted in equation (2), we foimd that a good correlation between the derived eq\ia- 
tion (1) and the data required that cx 2 ~ 0 for the data from capsules 122 and 124. This 
indicates, as will be discussed later on, that bvibble motion was not an important trans- 
port mechanism for these two capsules. 

Some of the capsule 124 tests were conducted using neon (Ne) , instead of helium 
(He), as the sweep gas. In this manner we were able to test this particular fuel speci- 
men at different power -temperature combinations. 

As reported in previous reports (refs. 4 and 5), the capsule used for our sweep gas 
tests contains a heat transfer gas gap between the fuel specimen and the capsule wall. 
Thus a capsule containing a particular gas gap width would have a set fuel power - 
temperature relation. Changing the sweep gas from helium to neon drastically changed 
the gap conductivity and allowed us to run tests on capsule 124 at: ! 

(1) The same temperature, but widely different power densities, or ^ 

(2) The same power density, but different temperatures. 

Release rate data for capsule 124 is presented in tables V(d) and VI of appendix A. 

This data was not used to obtain the constants for equation (2) but does serve as a good 

check on the temperature and fission rate density dependences. 

88 

A modified form of equation (2) was used to calculate Kr release rates into the 

-3 

neon sweep gas. Modification was necessary because the recoil constant, 3.65x10 in 

equation (2), foimd for helium, had to be altered to account for the increased relative 

stopping power of neon. Usii^ values from Evans (ref. 9) the increase in the value of the 

recoil constant is a factor of 3.29 when changing from 99. 8 percent He-0.2 percent A to 

90 percent No-0.2 percent A (neon gas composition). 

88 

For the Kr release into neon calculations, the modified form of equation (3) be- 
comes 
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Figure 3. - Comparison of measured and calculated Kr®® activities in neon sweep gas. 

j 

where there is no bubble motion term. ( 

88 

Calculated and measured Kr release rates into neon are shown in figure 3. The 
mean relative difference is 10.6 percent, and the' standard deviation of the relative dif- 
ference is 13.7 percent. 

The derivation of an equation to predict stable fission product gas releases from 
rectangular grains was carried out in reference 5. The result in terms of total release 
divided by total birth is 


00 



n=0 


9 


-Q./RT 


-Q^/RT 




(from eqs. (B5) and (B7)) and recoil release has been neglected. 

For eqviation (5) to be usefvil, values must be obtained for the diffusion constants 
Dfo, D^, and 02 - Values obtained for the diffusion terms in equation (1) contain the 
nonseparable surface area factor S^. 

If we assume, as was done in reference 5, that after approximately 6 percent burn- 

up we have attained complete interconnection of circumferential intergranular porosity 

2 4 

as in capsules 121 and 123 (fig. 1), S is approximately 200 centimeters , and D be- 
comes the following: 


18 800 


18 400 


D = 8.22xl0"^^f + 2.37xl0“^°e 


+ 1 . 00 x 10 '^® e 


( 6 ) 


For calcxilational purposes, equation (5) can be approximated by 


0.479tj^/^ for ?7 < 1 

(7) 

1 - 

2 2 

where 17 = Dtt t/ 4 y 2 . 



1,1 e-^ 

625tj 2401t7 t ] 


for Tj > 1 



Stable Fission Gas Release Correlation 

Equations relatii^ stable fission gas release from UN fuel pins to the release from 
single grains have been derived (appendix D) and have been correlated with available UN 
test data. 

The model used assumes fission gas release from individiaal grains to the grain 
boimdaries. If a grain is at the fuel surface, that fission gas is immediately released 
from the pin. Gas is only released from the fuel interior through interconnected porosity 
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which forms as the amount of intergranular gas increases. 

Two basic assumptions are that the intergranular gas is at a constant pressure, and 
that the fraction of porosity interconnected to the surface is proportional to porosity. 

The following equation (derived as eq. (D8)) relates total gas release from the pin 
R,p/B to total gas release from a grain (R,p/B)g: 



where P is the constant gas pressure and /3 is the factor which defines the rate at 
which porosity interconnects, is interconnected porosity /total porosity, and is found 
from the data analysis. (Rrp/B)^ is found from equation (7). 

The constant gas pressure P in equation (8) is assumed to equal the steady -state 
fuel- clad contact pressure which can be calculated from the following equation (derived 
as eq. (D(18) in appendix D); 


(nc-nf) Apk, 


4f 


A k. 

c 4c 


^4c^cP 


n. 


(9) 


The basic assumptions used in deriving equation (9) were that both the fuel and clad are 
creeping at the same constant rate, tangential and axial stresses in the fuel and clad 
are proportional to radial stress, and the fuel volumetric strain rate can be approxi- 
mated by the vinrestrained fuel swelling strain rate (calculated from Lietzke, ref. 10). 

The UN fuel pin data presented in appendix E were correlated with equation (8) . 
The.results of the correlation are shown in figure 4’ for those pins which had a measured 
gas release of over 1 percent. 

Calculated results for all the experiments are presented in table II. The average 
relative difference for all the data is 44.3, the standard deviation of the relative dif- 
ferences is 37 percent. _ . . 
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1 10 100 
Calculated gas release, percent of total formed 


Figure 4. - Comparison of measured and calculated stable gas release from UN (high 
gas release data only). 


For the data presented in figure 4, the average relative difference 



with a standard deviation of 0. 52. For this data the best value of ) 3 , fovind by a least- 
squares analysis, is 9.73. This indicates that the porosity is completely intercoimected 
when the total intergranular porosity reaches ~10 percent. 

Empirically it was also apparent that good agreement between equation (8) and the 
data required that bubble motion occurs only at the higher temperatures and burnup 
rates and that pore interconnection occurred only above some minimum temperature and 
burnup. At low temperature and burnup the gas release comes only from that layer of 
grains touchii^ the fuel surface. 
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TABLE II. - CALCULATED AND MEASURED STABLE GAS RELEASE DATA 


Experi- 

ment 

Fuel 

temper- 

ature, 

K 

Burnup, 

percent 

Burnup rate, 
percent 
10 000 hr 

Model^ 

Contact 

pressure, 

N/cm^ 

Rt/b. 

percent 

Meas- 

ured 

Calcu- 

lated'* 

1 

1144 

0.06 

1.54 

S 

944 

0.08 

0.031 

2 

1248 

.06 

1. 54 

S 

560 

.09 

.048 

3 

1443 

.057 

3.58 

B 

418 

3. 5 

5.4 

4 

1422 

.72 

4. 53 

B 

511 

5. 1 

5.8 

5 

1477 

.42 

3.85 

B 

390 

3.0 

3.6 

6 

1533 

.42 

3.85 

B 

327 

3. 3 

4. 4 

7 

1533 

.40 

3.67 

B 

324 

2. 5 

4.2 

8 

1429 

1.20 

7. 10 

B 

601 

4.8 

5.2 

9 

1582 

1.20 

7. 10 

B 

362 

11.7 

8.8 

10 

1443 

1.21 

7. 16 

B 

593 

4.2 

5. 5 

11 

1568 

.58 

6.82 

B 

398 

. 53 

6. 1 

12 

1658 

.63 

7. 41 

B 

314 

.88 

7. 5 

13 

1547 

.61 

7. 17 ■ 

B 

441 

.68 

5.3 

14 

1533 

2.26 

8.22 

B 

541 

9.9 

11. 1 

15 

1665 

.2.23 

8. 11 

B 

356 

29. 1 

18.2 

16 

1581 

2.34 

8. 51 

B 

00 

11. 1 

13.6 

17 

1554 

1.65 

5. 50 

B 

317 

17.2 

12.0 

18 

1637 

1.82 

6.01 

B 

268 

32.3 

17.3 

19 

1568 

1. 51 

5.03 

B 

283 

21.9 

12.6 

20 

1616 

.78 

7.09 

B 

285 

5.8 

8.8 

21 

1665 

.74 

6.73 

B 

240 

3.9 

10.3 

22 

1588 

.72 

6. 54 

B 

297 

4.0 

8.0 

23 

1348 

. 19 

2. 87 

S 

462 

. 1 

.037 

24 

1470 

. 18 

2.72 

S 

304 

. 1 

.095 

25 

1483 

. 18 

2.72 

S 

295 

. 1 

. 11 

26 

1403 

. 39 

2.-84 -- 

- s- 

392 

■ . 1 

.090 

27 

1470 

. 39 

2.84 

S 

314 

.2 

. 13 

28 

1458 

.40 

2.91 

S 

330 

.2 

. 14 

29 

1265 

.78 

2.35 

S 

690 

.05 

.081 

30 

1408 

.76 

2. 30 

S 

391 

. 1 

. 17 


^Interconnection, I; bubble motion + interconnection, B; surface- release, S. 
'interconnection-factor ^ = 8.75. - 




TABLE II. - Concluded. CALCULATED AND MEASURED STABLE 


GAS RELEASE DATA 


Experi- 

ment 

Fuel 

temper- 

ature, 

K 

Burnup, 

percent 

Burnup rate, 
percent 
10 000 hr 

Model^ 

Contact 

pressure, 

N/cm^ 

R,p/B, percent 

Meas- 

ured 

Calcu- 

lated*’ 

31 

1388 

0. 78 

2. 36 

S 

438 

0. 06 

0. 17 

32 

1240 

2. 54 

4.88 

S 

3760 

.05 

. 11 

33 

1350 

2. 54 

4. 88 

S 

2500 

.37 

. 17 

34 

1285 

2. 54 

4. 88 

S 

3120 

.09 

. 13 

35 

1475 

1. 18 

1. 18 

s 

975 

. 34 

.39 

36 

1460 

1.43 

1.36 

s 

1070 

.04 

. 38 

37 

1543 

. 70 

2. 20 

s 

991 

.09 

.29 

38 

1538 

.60 

1.89 

s 

963 

.03 

.29 

39 

1538 

.60 

1.89 

s 

963 

. 11 

.29 

40 

1538 

.60 

1. 89 

s 

963 

. 13 

.29 

41 

1480 

1.08 

1. 10 

s 

947 

. 14 

. 40 

42 

1460 

1. 08 

1. 10 

s 

996 

.09 

. 37 

43 

1558 

.07 

.65 

s 

698 

.05 

.28 

44 

1363 

1. 16 

1. 13 

s 

1330 

.01 

.24 

45 

1498 

1. 16 

1. 13 

s 

924 

.08 

. 44 

46 

1488 

1. 16 

1. 13 

s 

947 

.04 

.42 

47 

1228 

4. 58 

4. 78 

I 

4160 

1.36 

1. 50 

48 

1332 

4. 58 

4.78 

I 

2730 

5.96 

2. 37 

49 

1275 

4. 58 

4.78 

I 

3450 

1.78 

1.83 

50 

1473 

.78 

1.26 

s 

979 

.02 

.31 

51 

1478 

.78 

1.26 

s 

966 

.03 

. 31 

52 

1160 

2. 72 

2.63 

s 

4000 

.05 

.097 

53 

1300 

2. 72 

2.63 

s 

2170 

. 13 

. 18 

54 

1255 

2.72 

2.63 

s 

2620 

. 12 

. 15 

55 

1497 

1. 12 

1. 10 

s 

908 

.06 

.44 

56 

1443 

4. 17 

3.48 

B 

1640 

6.0 

9.9 

57 

1503 

4. 17 

3.48 

B 

1392 

10.3 

12.8 

58 

1528 

4. 17 

3. 48 

B 

1304 

11. 5 

14.2 

59 

1503 

4. 17 

3.48 

B 

1392 

12.6 

12.8 

60 

1393 

1. 78 

2. 16 

S 

4635 

. 14 

. 30 

61 

1598 

1. 45 

2. 50 

s 

676 

. 10 

. 59 

62 

1920 

1. 75 

3.01 

B 

506 

7. 10 

10.8 

63 

1255 

. 47 

. 58 

s 

1570 

.05 

. 14 

64 

1270 

.89 

1. 10 

s 

5200 

. 10 

. 15 

65 

1271 

.92 

1. 14 

s 

5220 

. 10 

. 15 


^Interconnection, I; bubble motion + interconnection, B; surface release, S. 
interconnection factor |3 = 8.75. 




Interconnection 

Surface release / r Bubble motion and 



Figure 5. - Regions of interconnection and bubble motion. 


The data indicate that; Bubble motion exists if the fuel temperature exceeds 1375 K, 

the bxirnup exceeds 0.4 percent, and the burnup rate exceeds 3 percent per 10 000 
0 

hours. Gas release is from the surface if there is no bubble motion or if either the fuel 
temperature is below 1225 K or the burnup is below 4 percent. 

The bubble migration, interconnection, and surface release regions are shown pic- 
torially on figure 5 . 


Discussion of Radioactive Gas Release 

The diffusional transport -release equations derived in appendix B and used to cor- 

go 

relate our online Kr release rate data only approximate the complex phenomena oc- 
curring inside fissioning fuel specimens. Not considered in the derivation are a broad 
rai^e of interactii^ mechanisms; bubble nucleation, bubble growth, bubble coalescence, 
and gas resolution from bubbles. Nor have we examined bubble motion to determine 
2 

These numbers are only approximations. They were obtained from the data in 
appendix E . 
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which of the bubble diffusion mechanisms (surface, volrnne, or gaseous) controls bubble 
velocity. 

We have assumed that fission gas transport in UN by bubble diffusion is primarily 
dependent on fuel temperature and temperature gradient (as in eq. (B6)), and we have 
also assumed that the diffusion coefficients describing each of the three diffusion mech- 
anisms (fission enhanced , atomic , and bubble) are additive in their effect on fission gas 
(as in eq. (B5)). 

In effect, we have reduced a very complicated problem to a relatively simple one 
(eq. (1)) which very adequately describes the radioactive gas release rates from our fuel 
specimens, as shown in figure 2. 

go O 

The Kr release rate data in figure 2 covers an activity range of from 3x10 to 

fi O 

4x10 disintegrations /cm -sec, and yet the average absolute relative difference 

A - A 

k, measured k, calculated 

A 

k, measured 

is only 13. 8 percent. This 13. 8 percent relative difference compares very favorably 
with the estimated imcertainties in activity measurements (±10 percent) and in the calcu- 
lated fission rate densities (±15 to 20 percent). 

Gas sample results for capsule 124 (table VII of appendix A) offer further verifica- 
tion of equation (1) which predicts that the release rate should be approximately pro- 
portional to For any two of the four isotopes measured (Kr®^“^, Kr®®, or 

Xe^®®) the relative activities in the sweep gas would be 

Activity j 

Activityg Yglxgj ^-Xgtjj 

where recoil has been neglected. Calculated and measured activity rates are compared 
in table III. 

go 

The best correlation obtained between our online Kr release rate data and equa- 
tion (1) required that the bubble motion term 

-Q^/R 




A 
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TABLE III. - COMPARISON OF CALCULATED AND MEASURED 


ISOTOPIC ACTIVITY RATIOS 


Isotope pair 

Calculated activity ratio, 
activity j^/activity2 if 
mechanism is - 

Measured activity ratio 

1 

2 

Capsule 124^ 

Capsule 123^ 

'“gf 

m 

r * f{\) 

Kr85m 

00 

0.487 


0.262 

0. 520 

0. 584 

KrSSm 


. 331 


.262 

. 349 

. 371 

Kr®® 

Kr^ 

1.48 


1.01 

1. 49 

1. 57 

CxeI85 

00 

1.91 


.717 

1. 45 

1.85 


^Average of 26 gas samples for capsule 124 (table VII) . 

'^Reported in ref. 5. 

c 135 

Neglects neutron absorption in Xe and iodine decay to xenon outside 


fuel. 


be set to zero for capsules 122 and 124. The postirradiation photomicrographs of cap- 
sules 122 and 124 (figs, 6 and 7) certainly do not show the large scale bubble formation 
and motion indicated in figure 1 for capsule 123. As indicated in table I, capsule 123 
was operated to 1512 K while capsules 122 and 124 were only operated as high as 1225 
and 1357 K, respectively. 


Discussion of Stable Gas Release 

The concept of fission gas release through interconnected grain boimdary porosity 
is not new. It has been used by many authors to explain gas release and swellii^ be- 
havior in nuclear fuels. Ainscorgh (ref. 8) shows how such intergranular porosity de- 
velops and interconnects in UO 2 fuel, 

Our use of an interconnection model is primarily based on the evidence of postirra- 
diation photomicrographs such as figures 1, 6, and 7. For one of our pins (capsule 122, 
fig. 6) there is little evidence of any intergranular porosity; for another (capsule 124, 
fig, 7) there appears to be some slight oriented porosity; and for another (capsule 123, 
fig. 1) there is extensive orientation and interconnection of intergranular porosity. 

The maximmn temperature at which each capsule operated increased from 1225 K 
for capsule 122 to 1357 K for capsule 124 to 1512 K for capsule 123 (table I). The higher 
the temperature of operation, the more intergranular porosity - at least for our fuel 
specimens which were operated to high burnups. 


17 

















Inside Outside 

radius radius 

Figure 6. - Section of fuel pin in capsule 122. X500. 


The actual derivation of the equations to predict stable gas release (as carried out 
in appendix D) contains a number of implicit and explicit assumptions. Oiir intent in the 
derivation was to include the important physical phenomena while obtaining relatively 
simple final relations for ease in interpretation and analysis. 

Chir basic assumptions are as follows: 

(1) Fission gas release is to the grain boundaries. This is an obvious assumption 

for the ease where intergranular porosity exists and the intergranular space acts a sink. 
We assume no mechanism operates to get the gas back into the matrix. Resolution is 
neglected. When there is little intergranular porosity and little gas release, the diffu- 
sion or nondiffusion of gas across a grain boimdary is umimportant. If gas could freely 
cross a boundary, then grain half -thickness in eq. (5)) would be replaced by the 

much larger dimension of half the fuel thickness and the calculated gas release would 
still be low. 

(2) The intergranular gas pressure is at a constant pressure P (related to the gas 
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Inside Outside 

radius radius 

(a) X500. (b) X250. 


Figure 7. - Section of fuel pin in capsule 124. 


volume through the perfect gas law) and is assumed equal to the fuel -clad contact pres- 
sime. 

(3) The fraction of intergranular porosity interconnected to the outside is directly 
porpiortional to the total intergranular porosity . Physically we would expect that the 
degree of interconnection would vary nonlinearly with porosity - starting out near zero 
until some critical porosity was reached - then increasing rapidly with additional gas 
volume - asymptotically approaching unity. A more complete description of this com- 
plex process, with various bubble sizes and velocities as well as competing gas trans- 
port processes would only becloud our analysis. 

(4) All the porosity is intergraniilar . There is no simple way of distinguishing be- 
tween inter and intragranular porosity as a time varient problem. 

(5) A reasonable calculation of fuel-clad contact pressure can be made by assuming 
steady -state and equal tangential creep rates for the fuel and clad, proportional stresses 
in both the fuel and clad, and a very simplified fuel swelling rate. 

The assumptions appear to work. As shown in figure 5 and table in, calculated and 
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measured gas releases agree very well - if bubble motion and interconnection are limi- 
ted to definite temperature, burnup, and burnup rate regions as in figime 5. The mini- 
mum bvirnup limits shown in figtire 6 do, in effect, add some nonlinearity to the porosity 
interconnection process, partially compensating for aforementioned assumption (3). 

There are some inconsistancies in the basic gas release data used in o\ir analyses 
and listed in table IX of appendix E. Comparisons of three of the tests (i. e. , 6, 12, 
and 21) are shown in table IV. 


TABLE IV. - INCONSISTENT GAS RELEASE DATA 


Test designa- 
tion 

Fuel tem- 
perature, 
K 

Time, 

hr 

Burnup, 

percent 

Measured 

Rt/b 

12 

1658 

850 

0.63 

0.009 

21 

1665 

1100 

.74 

. 04 

6 

1533 

1090 

. 42 

.03 

12 

1658 

850 

.63 

.009 

21 

1665 

1100 

.74 

.04 

6 

1533 

1090 

. 42 

.03 


In the first case (12 and 21), a factor of 5 difference in gas release is unjustified. 
In the second case (6 and 12), the much higher temperature has the lower gas release. 
In the third case (21 and 6), widely different temperatures result in similar releases. 

We have no explanation for the inconsistent data. It is mentioned only to demon- 
strate a common problem in attemptii^ to analyze fission gas release data. 


SUMMARY OF RESULTS 

A model has been developed to predict fission gas release from tmanium nitride 
(UN) fuel pins clad with various materials. When correlated with measured gas re- 
leases from pins primarily clad with Nb-lZr or PWC-11, operated to ~1650 K, the 
average deviation between calculated and measured gas release is 44 percent. The gas 
releases vary from 0.01 to 32 percent. If only the measured releases over 1 percent 
are considered, the average deviation is 26 percent. 

In the model, fission gas is released to the grain boundaries where the porosity 
eventually interconnects and links to the outside. Interconnection does not take place if 
the fuel is below 1225 K or below 4 percent burnup. Bubble motion, one of the mech- 
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anisms for intragranular transport, exists only if the fuel temperature exceeds 1375 K, 
and the burnup exceeds 0,4 percent and the burnup rate exceeds 3 percent per 10 000 
hovirs. Diffusion coefficients for the intragranular transport have been obtained from 
data measured in a sweep gas facility. 

Fission gas release rates from UN specimens (sweep gas experiment data) irra- 
diated to over 8 percent burnup at temperatures to 1789 K can be predicted to within 
14 percent (mean absolute relative difference) T^ere the release rate varies over three 
orders of magnitude. The standard deviation is 11 percent, 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 12, 1973, 

503-05. 
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APPENDIX A 


FISSION GAS RELEASE RATE DATA FOR UN FUEL SPEQMENS 

Fission gas activities in the sweep gas were determined from online scintillation de- 
tector counts or film physical samples of the sweep gas v^ich were analyzed by the Plum 
Brook Radiochemistry Section. Details of the procedures and the data for capsules 121, 
122, and 123 have been previously presented in reference 5. 

go 

All of the online and gas sample Kr release rate data used to obtain values for the 
unknown constants in equation (1) are presented in tables V and VI. 

Complete listings of the gas sample results for capsule 124, incomplete in time for 
inclusion in reference 5, are given in table VH. 


TABLE V. - KRYPTON-88 ACTIVITY IN HELIUM SWEEP GAS 


(a) Capsule 121 


Experi 

ment 

Fuel max- 
imum tem- 
perature , 

K 

Fission-rate 
density , 
fissions 

cm^-sec 

Fuel 

burnup, 

percent 

Type of 
data^ 

K®® activity, 
disintegrations 

cm^-sec 

Measured 

Calculated*^ 

1 

1550 

5. 77X10^“^ 

1. 50 

0 

554x10^ 

615x10^ 

2 

1550 

5. 77 

1. 83 

0 

720 

664 

3 

1550 

5. 77 

1.99 

o 

720 

687 

4 

1550 

5.77 

2.23 

0 

860 

723 

5 

1787 

8.28 

3. 46 

o 

2810 

2710 

6 

1787 

8. 28 

5.06 

o 

3630 

3420 

7 

1765 , 

8. 10 

5.06 

s 

2560 

3160 

8 

1778 

8. 18 

6.03 

o 

3970 

3710 

9 

1778 

8. 18 

6.59 

0 

4240 

3950 

10 

1765 

8. 10 

6.65 

s 

3690 

3810 

11 

1778 

8. 18 

6.87 

o 

4380 

4070 

12 

1218 

3. 34 

8.07 

- o 

215 

228 

13 

1215 

3. 34 

8. 30 

s 

183 

230 

14 

1218 

3. 34 

8.77 

o 

263 

243 

15 

1218 

3. 34 

8. 84 

0 

288 

244 

16 

1218 

3. 34 

8. 93 

o 

285 

246 


^Sample, S; online, O. 

^Value calculated using eq. (3). 
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TABLE V, - Continued. KRYPTON-88 ACTIVITY IN HELIUM SWEEP GAS 


(b) Capsule 122 


Experi- 

ment 

Fuel max- 
imum tem- 
perature, 
K 

Fission-rate 

density, 

fissions 

Q 

cm -sec 

Fuel 
burnup , 
percent 

Type of 
data^ 

activity, 

disintegrations 

cm^-sec 

Measured 

Calculated'^ 

1 

1223 

3.64X10^^ 

2.64 

S 

73.5 

65. 2x10^ 

2 

800 

1. 42 

2. 95 

O 

14. 3 

17. 4 

3 

996 

2. 33 

3.02 

O 

29.3 

34. 9 

4 

1219 

3. 59 

3.09 

o 

68. 7 

69.2 

5 

598 

.71 

4. 17 

o 

7.00 

8.61 

6 

805 1 

1. 44 

4.22 

o 

16. 5 

21. 8 

7 

1002 

2. 37 , 

4.31 

o 

34. 3 

43. 8 

8 

1000 ^ 

2. 33 : 

4.31 

s 

38. 5 

42.8 

9 

1223 

.3.65 

4.40 

o 

86.8 

87.2 

10 

1225 

3.64 

5.38 

s 

104 

99. 1 

11 

594 

.692 

5.40 

o 

9. 30 

9.82 

12 

804 

1. 44 

5. 44 

o 

26.9 

25.6 

13 

1002 

2. 37 

5. 54 

0 

46. 6 

51.6 

14 

1223 

3.65 

5.64 

o 

102 

103 

15 

598 

. 71 

5.73 

o 

9.60 

10.6 

16 

805 

1. 44 

5.81 

o 

25. 2 

26. 8 

17 

800 

1. 44 

5.81 

s 

25.7 

26.8 

18 

1000 

2. 33 

5. 93 

s 

57. 2 

52.7 

19 

1000 

2. 35 

5. 93 

o 

47. 3 

53. 4 

20 

1220 

3.60 

6.03 

s 

101 

105 

21 

598 

.71 

6. 16 

o 

8.00 

11. 1 

22 

805 

1. 45 

6.27 

s 

26. 1 

28. 5 

23 

1227 

3.64 

6.59 

- S ■ 

118 

114 


^Sample, S; online, O. 

'’value calculated using eq. (3). 
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TABLE V. - Continued. KRYPTON-88 ACTIVITY IN HELIUM SWEEP GAS 

(c) Capsule 123 


Fuel 

burnup 

percent 


Fuel max- 
imum tem- 
perature, 

Fission-rate 
density , 
fissions 

K 

cm^-sec 

585 

0. 505x10^^ 

^95 

1. 10 

795 

1. 10 

990 

1.83 

1228 

2. 79 

1387 

3.60 

1542 

4. 44 

580 

.496 

580 

.496 

803 

1.09 

803 

1.09 

995 

1.77 

995 

1. 77 

591 

. 505 

591 

. 505 

591 

.505 

804 

1.09 

804 

1.09 

1001 

1.88 

991 

1.84 

991 

1.84 

1222 

2.77 

1222 

2. 77 

1380 

3. 56 

1380 

3. 56 

1541 

4. 44 

1541 

4. 44 

585 

. 505 

995 

1.77 

1223 

2. 77 

1390 

3.62 


I 


Type of 
data^ 

K®® activity, 
disintegrations 


cm^ 

-sec 

Measured 

Calculated*’ 

0 

3. 85X10^ 

3. 33X10^ 

0 

9.65 

9. 07 

s 

10.7 

9.07 

0 

23.7 

20. 1 

0 

61. 0 

73. 5 

o 

152 

197 

o 

480 

467 

o 

5.35 

4. 87 

o 

5. 35 

4. 87 

s 

16. 4 

13. 4 

0 i 

14.8 

13.4 

o 

35.0 

28.9 

s 

44. 1 

28.9 

0 

4.30 

5. 57 

s 

6.68 

5. 57 

0 

5.00 

5. 57 

o 

17.0 

15. 1 

s 

22.0 

15. 1 

0 

45.0 

35.3 

o 

42.0 

35.9 

s 

46.8 

36.6 

0 

130 

129 

s 

108 

129 

0 

350 

343 

o 

350 

343 

0 

950 

836 ' 

s 

939 

836 

0 

6.40 

6.82 

o 

45.0 

38.8 

0 

175 

142 

0 

470 

393 


Sample, S; online, O. 

Value calculated using eq. (3) 



















Experi 

ment 


1 

2 

1 ^ 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 


Fuel max- 
imum tem- 
perature, 
K 

Fission-rate 

density, 

fissions 

cm^-sec 

Fuel 

burnup, 

percent 

Type of 
data^ 

K®® activity, 
disintegrations 

3 

cm -sec 

Measured 

Calculated*^ 

595 

1.01X10^'^ 

1. 53 

o 

8. 69x10^ 

8. 57x10^ 

594 

1. 10 

1.54 

S 

9.30 

9.67 

1002 

3.81 

1.84 

0 

53.2 

56.3 

1230 

6. 18 

2.30 ! 

o 

123 

130 

1230 

6. 18 

2.38 ' 

o 

148 

132 

808 

2.25 

2.58 

o 

29.2 

31.3 

810 

2.35 

2. 58 

s 

29.2 

33. 3 

425 

.357 

2.59 

o 

3.64 

2.99 

570 

.904 

2.66 

o 

11. 4 

9.45 

480 

. 535 

2. 66 

o 

4. 79 

4.91 

594 

1. 10 

2.66 

s 

14. 5 

12. 3 

836 

2.45 

2.85 

o 

42.9 

37.0 

1306 

7.28 

3.33 

o 

212 

208 

1296 

7. 13 

3.53 

o 

205 

207 

1294 

7. 12 

3. 53 

s 

1 

193 

206 

1002 

3.81 

3. 55 

-o 

69. 5 

77. 7 

632 

1.21 

3.66 

o 

30.2 

16.4 

626 • 

1. 16 

3.66 

s 

30.2 

15.4 

900 

2.96 

3.89 

. o 

53.9 

57.3 

1316 

7.45 

4.48 

o 

252 

258 

1323 

7. 57 

4. 49 

s 

241 

266 

592 

1.09 

4. 52 

s 

20. 9 

16. 3 

594 

1.01 

4.52 

o 

17. 8 

14.7 

594 

1.01 

4. 52 

o 

20. 6 

14.7 

1003 

3.82 

4.87 

o 

97. 1 

94. 5 

1229 

6. 17 

5. 39 

o 

229 

210 

1232 

6.33 

5.40 

s 

202 

218 

439 

.389 

5. 90 

o 

4. 22 

5.38 

.,596 - - 

1.01- - -- 

5.98 

o 

36. T 

17.7 

595 

1. 10 

5.99 

s 

38.6 

• 19.9 

810 

2.26 

6.01 

o 

77. 8 

52. 1 

810 

2. 26 

6. 16 

o 

64.2 

53.0' 

1016 

4. 13 

6.48 

s 

109 

128 

1012 

3. 93 . . 

6.48 

o 

121 

119 


6. 14 

7.53 

o 

229 

263 


7. 50 

7.75 

o 

444 

381 

1318 

7. 51 

7.76 

s 

454 

382 


















Experi- 

ment 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 


Fuel max- 
imum tem- 
perature, 

K 

Fission-rate 

density, 

fissions 

cm^-sec 

Fuel 

1 burnup, 
percent 

Type of 
data^ 

K®® activity, 
disintegrations 

cm®-sec 

Measured 

Calculated' 

849 

l.OlxlO^^ 

1.59 

0 

12. 0x10® 

11.3 

851 

1.07 

1.60 

S 

10.7 

12. 1 

999 

1. 55 

1.64 

o 

19.9 

18.8 

1219 

2. 50 

2.50 

0 

50. 1 

44.0 

1219 

2.49 

2. 51 

s 

54.4 

43.8 

429 

.21 

2.60 

o 

3.64 

3.20 

570 

. 39 

2.61 

0 

5.33 

5. 34 

595 

.43 

2.61 

s 

5.02 

5.90 

841 

1.00 

2.72 

0 

16.0 

14. 1 

1219 

2. 50 

2.90 

0 

74.2 

47. 3 

1358 

3. 19 

3.08 

0 

no 

80. 1 

1363 

3. 15 

3. 10 

s 

113 

79.9 

629 

. 48 

3. 59 

0 

7. 55 

7.60 

631 

. 50 

3.60 

s 

7.39 

8. 04 

895 

1. 18 

3.72 

0 

19. 5 

20. 2 

1303 

2. 90 

4.06 

0 

92.8 

76.0 

1308 

2.88 

4.06 

s 

105 

75. 8 

833 

.99 

4. 57 

o 

21. 5 

18.8 

835 

1. 06 

4.58 

s 

22. 4 

20.6 

999 

1. 59 

4.65 

s 

41. 1 

33. 1 

1000 

1. 55 

4.65 

o 

36.9 

32.0 

1225 

2. 53 

5.53 

0 

75. 0 

71. 2 

1221 

2. 50 

5. 54 

s 

81. 8 

70. 0 


















TABLE Vn. - FISSION GAS ACTIVITY FOR CAPSULE 124 (GAS SAMPLES) 

(a) In helium sweep gas 


Fuel 

burnup 

percen 


Fuel max- 
imum tem- 
perature, 

Fission-rate 

density, 

fissions 

K 

cm^-sec 

806 

2.33X10^^ 

803 

2.29 

554 

.90 

1008 

3.99 

594 

1.10 

810 

2.35 

594 

1.10 

1294 

7.12 

626 

1.16 

1323 

7.57 

592 

1.09 

1232 

6.33 

595 

1.10 

1016 

4.13 

1318 

7.51 


Isotope activity, disintegrations/c 

m -sec 

85m 

Kr 


Kr®8 

Xe^®® 

7.96x10^ 

18.2x10® 

25.2x10® 

22. 2X10^ 

8.16 

16.5 

24.6 

21.4 

2.07 

3.35 

6.02 

6.33 

19.6 

39.3 

56.5 

48.3 

3.23 

5.93 

9.30 

8.96 

10.0 

19.0 

29.2 

28.5 

5.02 

8.95 

14.5 

14.7 

69.1 

129 

193 

195 

10.8 

20.0 

30.2 

30.4 

87.9 

162 

241 

238 

7.62 

14.2 

20.9 

20.3 

70.3 

138 

202 

195 

13.5 

28.3 

38.6 

38.6 

38.1 

71.6 

109 

103 

173 

278 

454 

436 


(b) In neon sweep gas 


554 

0.35X10^^ 

0.84 

0.886X10® 

1.88X10® 

2.73X10® 

2.44X10® 

999 

1.59 

1.11 

7.96 

16.4 

23.3 

20.4 

851 

1.07 

1.60 

3.85 

8.59 

10.7 

10.8 

1219 

2.49 

2.51 

19.8 

36.4 

54.4 

57.6 

595 

.43 

2.61 

1.65 

3.45 

5.02 

5.02 

1363 

3.15 

3.10 - 

.40.5 - - 

75.5 

113 

128^ 

631 

.50 

3.60 

2.62 

5.47 

7.39 

7.85 

1308 

2.88 

4.06 

35.5 

67.9 

105 

107 

835 

1.06 

4.58 

7.85 

15.0 

22.4 

23.7 

999 

1.59 

4.65 

14.9 

27.5 

41.1 

43.2 

1221 

2.50 

5.54 

29.0 

58.5 

81.8 

86.9, . 























APPENDIX B 


FUEL SPEQMEN FISSION GAS RELEASE RATE EQUATION 

One -dimensional diffusion of a radioactive fission product gas through a fuel grain 
in a rectangular coordinate system can be described by the following equation: 

2 

^=D2-S.-xc+fY (Bl) 

where the diffusion coefficient D is assumed to be constant, and the fission-rate den- 
sity f is independent of position in the fuel. 

In steady state with 9c /9t = 0 and with botmdary conditions 9c /9y = 0 at y = 0, 
and c = 0 at y = y^, the solution of equation (Bl) is 


c(y)=Sf 

X 



The fission gas release rate from the grain siurface can be foimd by adding a direct 
recoil component (ref. 5) to the dtffusional release, so that 


R = - DS —I 

S 4 dy 


(B3) 


Differentiatir^ equation (B2) and applying the result to equation (B3) gives 


O 




which in turn reduces to 



PYSLf 



4 


(B4) 



yi»i 



(which holds for all cases of practical interest). 

If it is now assumed that there are three additive diffusion transport mechanisms 
(fission-enhanced diffusion, atomic diffusion, and pore migration) and furthermore that 
all pores have a constant velocity and diameter, D can be replaced by 

D =D,f + De ° + V d„ (B5) 

fo o p p ' ' 


Nichols in reference 2 shows that the velocity of a pore whose motion is controlled 
by surface diffusion can be given by (in his notation) 



where 

D surface diffusion coefficient 

V number of rate -controlling species per unit area 

s 

molecular volume 

* 

Q„ surface heat of transport 

s 

k Boltzmann constant 

T absolute temperature 

r pore radius _ _ _ _ - _ 

K matrix thermal conductivity 

K' pore thermal conductivity 

VT^ macroscopic temperature gradient 

For our applications to cylindrical geometries we can replace equation (B6) first by 
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V d 
P P 




-Qg/RT 



where is a constant, and then substitute dT/dr at the outer fuel radius for dT/dy 
where 


dr 


^2 


2. 75x10"^ ^fr 2 


2K 


1 



(fuel power in W/cu cm = 2. 75x10 


"^^f) to get 


^ ^ pYSLf ^ fYS 
e" 4 ,1/2 


-Q„/RT -Q,/RT 

IDfof + Doe ° taje = 


a/2 




(B7) 


where 


“2 = 


2.75X10 


-11 



The recoil constant p is approximately proportional to gas density, so it may be 
approximated by Og/T, where «g is a constant. As seen in reference 5, the experi- 
mental gas release rate data shows Rpjj^ increasing linearly with burnup, so the final 
form of equation (B7) becomes 



= (1 + a^b) 



f -Q_/RT a„f -Q„/RT^ 

Df+De ° “ 

1/2 

a,YSLf 
^ + fYS 



4T 

X 

1 



V / 



(B8) 
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where b > 1 . 5 percent and s now refers to the total surface area through which gas is 
being released to the outside. 

In our analysis we used the following simple relation between thermal conduc- 
tivity Kj and temperature T: 

Kj(W/cm-K) = 0.1227 + 1X10~‘^T 
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APPENDIX C 


NONLINEAR LEAST -SQUARES CORRELATION 

The method of least sqiiares can be applied to an equation nonlinear in terms of un- 
known constants. The procedure is presented in Scarborough (ref. 11) and is also out- 
lined here. 

(1) Assume you have an equation y = f(x,aj,a 2 ,a 2 ) which is nonlinear in the un- 
known constants a^, ag, and a^, and a set of measured data points 

(x. ,y.,), (x„,y„), . . . , (x_,y ). Here y and x are general dependent and inde- 
pendent variables. 

(2) Choose initial values for a^, ag, and a^, such as a^Q, agQ, and agQ where we 
would have 


^1 ' ^10 h 
^2 "" ^20 ^2 


^3 " ^30 ®3 

when 6j, 6g, and 6g are unknown corrections. We are going to iteratively converge 
on values for a^, ag, and which minimize in a least-squares sense the differences 
between calculated and measured values of y. 

(3) Expand the equation y = f(x,aj^,ag,ag) about the initial values a^Q, agQ, and 
a^Q and retain only the first derivative terms so that, for each data set i, 

/df.\ /0f.\ 



(4) The "residual" equations of the following form: 

V. = y calc i - y. 

j j •'i 


are 





i = 1, N 
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For the method of least squares we want to minimize 


N 

g = ^ (Ci: 

i=l 

To do so, we obtain equations for 0g/36j, dg/dd^, and 9g/36g and set each equal to 
zero, and obtain three equations which can be solved for the uiicnown corrections 6j, 

62 , and 6g. 

(5) Next, return to step (2) and obtain the following better values for the constants 
aj, ag, and a^: 

^11 = ^10 ^1 
^21 " ^20 ®2 
^31 "" ^30 ®3 

and repeat the entire process as often as needed until there are no significant changes in 
aj, ag, or ag 

The followii^ equation (derived as eq. (B8)) was correlated with our Kr release 
rate data using this method: 


Rpin = (1 - «4^) 


4. f n ^ “2^ -Q./Rt\ 


1/2 


dgYSLf 

4T 


+ fYS 


Dr„f + D e 
fo o 


k|t2 


(C2) 


Values were obtained for the seven \inknown constants a^, a^, D^, Q^, cng, 

and Qg. Some specific steps are presented next. 

: : : To prevent computational problems,- equation (G2) was scaled and-rewritten as 
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r 




.,3 „ f,2 23 




1/2 




(C3) 


where the unknown constants are a^ to a,^. f, and ^ have been defined as 




-Xt, 


10 ^ 10 ^ 


f 


f’ “ 


10 


14 


A T 
S~ = 


10 ‘ 


oo o ^ 

Measured Kr activities vary over three orders of magnitude (~3xl0 tO'«4xl0 

Q 

disintegrations /cm -sec), so in an effort to keep the higher activity data points from out- 
weighing the lower activity points in the data fitting procedure, we minimized equations 
of the form 



instead of 


g = 



where 


^i '®^i, calc Vila ted 



measured 
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In expanded form g is 



and the seven equations in the seven imknown corrections 6^ which are found by settii^ 
0g/06^ = 0,i = l, . . . , 7, are of the form 

n 

■I 

i=l 

Application of this fitting procedure to the data of table V of appendix A resulted in 
the equation 
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appendix d 


STABLE GAS RELEASE FROM FUEL PIN 

The extention of our previously developed fission gas release equations (appendix B) 
to describe the release of stable fission product isotopes from reactor -type fuel pins is 
accomplished in two steps. First eqiiations relating gas release to intergranular poros- 
ity and intergranular gas pressure are developed and then the intergranular gas pres- 
svire is related to fuel and clad strer^th properties. 


Porosity and Stable Gas Release 

If it is assumed that the intergranular fission gas is at a constant pressiire P, the 
rate of gas volume increase is given by the voliune balance equation 


dV 


£ = rn 

dt S g p p ® p 


(Dl) 


. Then y is defined as the fraction of intergranular gas volume connected to the fuel 
surface and so the release rate from interconnected porosity is represented by 


R. = (R N - R)y 
ip g g 


so equation (Dl) becomes 


dV 


-S. = R N - (R N - Rjy 
dt g g p ' g g 


CT _r ^ 

P ® P 


(D2) 


Assume that y is proportional to the intergranular porosity 

V 

r = (3 (D3) 

Vg.V, 

where the coefficient /3 governs the rapidity at which interconnection occurs. As an 
approximation 
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and equation (D2) then becomes 


y 



^ = (R N - R ) - (R N - R ^ — 

dt g e ® P s g p 


or 


(R N “ R ) T^rp ic’v 

dv + ^ 0V dt = (R N - R ) ^ dt 

g Vf P^ ggsp 


The solution to this linear first-order equation is (ref. 12) 


V„ exp 
g ^ 



g i 


R ) dt 


PV^ 





x<exp 



R ) EL§. dt 


PV, 


f J 


dt> I + constant 


(D4) 


The integration 


- /* (R-N - R ) Rl^dt ' - - 

y ' g g s' pvj 

in equation (D4) can be performed if we have both R and R„ as functions of t. In our 

& ® 

previous report (ref 5) we showed that the diffusional release from a rectangular grain 
(neglecting svirface recoil) is 
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Rg = fYy2S 



exp< -D 


(2n + 1) j7_ 

2 y2 


-I 2^ 

t 


(2n + 1)2 


2 2 3 

For values of Dtt t/4y2 < 1 equation (D5) can be approximated by 


R =0.718 fYy2S 


fD^t 

4yi 


L/2 


(D5) 


or 


R = 0.359 fYV 

o o 



t = dt^/2 


if we also assume that 


R = FR N„ = FCN t 
s g g g 


1/2 


equation (D4) becomes 


V„ exp I y* (1 - t ^/2 KT^dtl = 


S PV, 


f 

f 


(1 - F)dN 

g p 


X 


jexp 


/ 


(1 - F)dN t^/2 KT^dt 


g 


PV, 


dt > 1 + constant 


^This corresponds to an instantaneous release to birth rate of less than 70 percent. 
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which yields 


exp 

g 


2 KTg(l - DdN^ ^3/3 

3 PV, 



(1 - F)dN t'/2 KT 
e p 


X J exp 


2KTg(l - F)dN^ ^3^2 


3PV, 



+ constant 


If we remove V,//3 from the integral on the right, it becomes of the form 

/ 6 ^6 

u'e dt with solution l/Se®'^. Thus the solution becomes 


V„ exp 
g 


„ KT^(1 - F)dN t^/2 


PV. 


f 

= — exp 

/3 


2 KT^(1 - F)dN t 


3/2 


PV. 


+ constant 


or 


V =-L 



r~ n 


/ 1 + h exp 

„ KT0(] - F)dN t^/2 
^ g 

K 

3 PVj 



- 

J 


Then 


The constant h can be determined from the initial condition V„ = V at t = 0. 

g go 


Y=-L 


' \ 

1 4- 1 — 1 1 

2KTj3(l - F)dN t^/^ 

o 

-N 

y 1 V,- 

3PVj 


\ / 

_ — ^ 

> 


(D6) 


The fission gas release rate from a fuel pin can now be calculated as follows; 
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/3V, 


+ Ro = (R N - R Jy + = (R N - R J — ^ + R 

pin ip s ' g g s ' g g s' „ 


Using equation (D6) and R = FN R we get 

& s 


Vg^/3 


2KT^(1 - F)dN t^/^ 

^ Q 


3PVr 


+ FR N (D7) 

g g ' 


1 /2 

If we rearrange equation (D7) and replace R by dt ' , we can integrate to get the 

& 

total gas release from a pin 




f) 

r 

exp 

2 KTj3(l - F)dN t^'^2 

"S 

dtl 

/ 

3 PVj 

/o 

«... 

- 



nin = -dN t^/^ 
T,pin 3 g 3 


f 


KT/3 


' 1 - exp 


2 KT^d - F)dNg t3/2|l 


3 PV, 




(D8) 


Equation (D8) can be put into terms of total release to total birth by dividing by 
fYVjt to give 


2 3/2 2PVf— ^-11 

P ^dN H Vf 

3 g \ f 


B 


fYVjt 


1 - exp 


3KT^Vjt 

Since the total release from a grain is given by 


2 KT|3(1 - F)dN t 


3/2 


PV. 


(D9) 
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dt = ^ 
3 


and the total birth in a grain is fYV t, equation (D9) becomes 

o 



Xjl-Ar 


2P - 1 


► ■ 


/HtV }\ 

1 - exp 

KT0(1 - F)fYt 

V 

) 3KT/3fYt 

■ ^ Iv 




' 'si 



For a cylindrical fuel pin with an outside radius rg and an inside radius 
F = 2y2/(r2 ~ final form of the release equation becomes 



Jl -exp 


KTjSfYt 



> 


(DIO) 


In doing the calculations with equation (DIO), the quotient V /V, is considered to 

&0 ^ 

be the initial fuel porosity and (Rrr,/B)_ is calculated from the integrated form of equa- 

^ 1/2 

tion (D5), not the approximate form R = dt ' . 

o 


■ Fuel-Clad Contact Pressure - - - - - 

The tai^ential creep rate (fuel or clad) is given by (ref. 13) 

= (Dll) 

e 

Assume that the steady -state equivalent strain rate and equivalent stress are re- 
lated by 
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e = Ao” 
e e 


where A is a function of temperature. Equation (Dll) becomes 




(D12) 


The equivalent stress can be written as 


cr = 
e 


fT2 


(-r - '^ef * * ("z - ‘’r)‘ 


nV2 


and if we assume the following proportional stresses: 


= k, cr 
9 1 r 




(D13) 

(D14) 


then becomes 


e ^ 1 r 


V~2 


(l-kjf.(k,-k2)2.(k2-l)2j 


1/2 


= k„ cr 
3 r 


(D15) 


where kg is a constant. Combining equations (D12) and (D15) yields 




(D16) 


where 


k4 = lk^"-l)(acj-k2-l) 


If we now equate fuel and clad tangential creep rates 


^6c = ^9f 
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where v is the fuel bulk strain rate, and substitute for using equation (D16), we get 



rc 




rc 


Afk4f 


‘^rf 


(nj-1) 

^rf 



(D17) 


At the interface we equate the radial stresses and fuel -clad contact pressure 
= -P, and equation (D17), after rearranging, becomes 


(Hc-nf) Ajk^j j 


A„k.^ 

c 4c 


n. 


^^4c^cP 


(D18) 


To use equation (D 18) , values for k^ and v have to be found. 

The constant k^, defined by equation (D13), can be approximated using the equations 
for stress in thick -walled cylinders (ref. 14) 


cr 

r 



- p„) 


(D19) 



C^P 


+ /^f(P. - P ) 
o I r/ ^ ° 


(D20) 


where c is the external radius and a is the internal radius. 

For the fuel, equations (D19) and (D20) are evaluated at the fuel outside radius 
with P^ = 0, giving 



and for the clad, equations (D 19) and (D 20) are evaluated at the clad inside radius with 
P^ = 0, yielding 



(D22) 
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Then k^, defined by equation (D14), can be approximated by an elastic calculation 
assuming that the fuel and the clad do not slip. 

In the axial direction 




and since 


a a , 
zc zf 




E Ef 
c f 


The forces creating and are equal and opposite , so 


•^'7 


S E " S„Ej E ^ ^ 

c c f f c f 


where and are the cross-sectional areas of the fuel and clad and 3^ is the 
force. 

Evaluating equation (D13) at the fuel-clad interface where 



yields 


(D23) 
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- 4 




IS E 
c c 


SfEf 


and since a = and a , = ^„/Sf, we obtain from equation (D24) 

ZC Z C ZI Z X 



Then can be calculated from kj and k 2 , since 

k4 = lk<'‘-‘>(2kj-k2-l) 


(D24) 


(D25) 


(D26) 


and 


k 


3 



An approximation to the fuel bulk strain rate' v can be obtained by considering the 
fuel swelling to be unrestrained. Lietzke (ref. 10) obtained integral equations for un- 
restrained swelling and has presented the results graphically with vt plotted as a 
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4 

function of a dimensionless parameter II for various values of the fuel creep rate 
stress exponent n^. 

Two simple empirical approximations to the integral equation result for vt are 


VT = 6.5 


(n,+ s) 


0.285jj^“f 


(nf+3)/n. 


r 

n£^3, n> 0.1, or 
n, > 4, n > 0.04, or 

< 

nj > 5, n> 0.02, or 
n^ > 6 


(D27) 


and 


VT = 





^(nj+l)/nj 


n^ < 3 , n < 0 . 1 , or 
nj < 4, n < 0.04, or 
nj < 5, n <: 0.02, or 
n^ < 6 


(D28) 


Comparisons of the simplified relations with the integral solution are shown in 
table Vm. 


l/(nj+l) ,Hj/(nj+l) 

( 3CyKTNB | 

\ "f / 

where in Lietzke's notation t is test time, B is fractional bxrrnup, N is fissionable 
atom density, and CY is fraction of gaseous fission products. 


n = 


3At 


2(uj + 1) 
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TABLE VIII. - UNRESTRAINED FUEL SWELLING 


[Comparison of approximate and integral solution for vt. ] 


"f 

TI 

0. 01 

0. 02 

0.04 

0.06 

0.08 

0. 10 

0. 12 


VT 

integral 

_ ^ _ 

0.01 

0. 027 

0.050 

0.080 

0. 117 

0. 160 

n 

VT 

approx. 

— 

.014 

.035 

.059 

.087 

. 108 

. 156 

4 

VT 

integral 

— 

.016 

.044 

.083 

. 136 

. 201 

. 28 

1 

VT 

approx. 

.... 

.019 

.040 

. 082 

. 136 

.201 

.28 

4 

VT 

integral 

0.01 

.023 

.064 

. 122 

.205 

.295 

>.33 

[ 

VT 

approx. 

.01 

.022 

.068 

. 130 

.207 

. 295 

>.33 

4 

VT 

integral 

.012 

.031 

.087 

. 175 

.295 

>. 33 


1 

vr 

approx. 

.012 

.034 

.097 

. 179 

.275 

>. 33 


4 

VT 

integral 

.015 

.039 

. 115 

.23 

>.33 



1 

VT 

approx. 

.017 

.047 

. 126 

.23 







APPENDIX E 


STABLE GAS RELEASE DATA FROM UN FUEL HNS 

All of the stable gas release data used to correlate with the equations developed in 
appendix D are presented in this section. The data comes from various sources (refs. 
15 to 19) and describes gas releases from pins clad with T-111, Nb-lZr, PWC-11, and 
W-26Re irradiated at temperatvires ranging from 1100 to 1900 K. 

Grain sizes for most of the pins listed in table DC were estimated from photomicro- 
graphs of fuel in similar capsules. 

Materials properties for the various clad materials along with the values used for 
UN are given in table X. 
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TABLE IX. -STABLE GAS RELEASE DATA 


E\ih-i 1 - 
mcnt 

Test 

description 

Fuel 

code® 

Clad 

code*^ 

Fuel 

temper- 

ature, 

K 

CUd 

temper- 

ature, 

K 

Time, 

hr 

Gram 

size, 

cm 

Measured 

R/B 

Clad 

outside 

radius, 

cm 

Clad 

inside 

radius, 

cm 

Fuel 

inside 

radius, 

cm 

Imtial 

porosity 

Fission-rate 

density, 

fissions 

cm^-sec 

Burnup, 

percent 

1 

PW26-200 42 


4 

U44 

1103 

390 

0.00330 

0.0008 

0.3985 

0.3100 

0 

0.096 

0 16Q0E^14 

0.060 

2 

PW26-200 43 





1248 

1214 

390 



.0009 







.090 

. 1430E-fl4 

060 

3 

PW26-201 40 





1443 

1353 

1 590 



.0350 







092 

.3530E4-14 

570 

4 

PW26-201 41 





1422 

1311 

1 590 



0510 







.097 

.4470E4-14 

720 

5 

PW26-210 108 





1477 

1380 

1 090 



.0300 

3 

925 

.3 

335 



.039 

.4040E-fl4 

.420 

6 

PW26-210 109 





1533 

1436 

1 090 



.0330 







.039 

.4070Efl4 

420 

7 

PW26-210 no 





1533 

1436 

1 090 



.0250 







038 

.3620E-fl4 

.400 

6 

PW26-220 118 





1429 

1297 

1 690 



.0480 







.059 

.5600E-fl4 

1.200 

9 

PW26-220 119 





1582 

1450 

1 690 



.1170 







.057 

.5560E+14 

1.200 

10 

PW26-220 120 





1443 

1297 

1 690 



.0420 







.060 

.5640E-fl4 

1.210 

11 

PW26-230 123 





1568 

1408 

850 



0053 

3 

960 

.3 

75 



.056 

.6620E^14 

.580 

12 

PW26-230 124 





1658 

1492 

850 



.0088 







.040 

.7090E+14 

.630 

13 

PW26-230 125 





1547 

1380 

850 



0068 







043 

.6910E-fl4 

.610 

14 

PW26-231 126 





1533 

1353 

2 750 



0990 







.040 

.7820E+14 

2.260 

15 

PW26-231 127 





1665 

1492 

2 750 



.2910 







039 

.7710E*14 

2.230 

16 

PW26-231 128 





1581 

1394 

2 750 



.1110 







038 

.8070E+14 

2 340 

17 

PW26-240 132 





1554 

1408 

3 000 



.1720 



3125 



.050 

S820E-fl4 

1.650 

18 

PW26-240 133 





1637 

1478 

3 000 



.3230 







.051 

.6550E4-14 

1.820 

19 

PW26-240 134 





1568 

1436 

3 000 



2190 







056 

.5460E-fl4 

1 510 

20 

PW26-241 135 





1616 

1436 

1 100 



0560 







051 

.7490E-fl4 

.780 

21 

PW26-241 136 





1665 

1492 

1 100 



.0390 







.053 

.7270E+14 

740 

22 

PW26-241 137 





1588 

1422 

I 100 



.0400 







.057 

.6980E+14 

.720 

23 

PW600 TOP 





1348 

1313 

662 



.0010 

3175 

.2415 



.022 

.1660E4-14 

.190 

24 

PW600 MID 





1470 

1438 

662 










.039 

1960E4-14 

.180 

25 

PW600 BOT 





1483 

1448 

662 










.045 

.1930E^14 

.180 

26 

PW601 TOP 





1403 

1368 

1 372 










.033 

niOE+14 

.390 

27 

PW601 MID 





1470 

1438 

1 372 



.0020 







.038 

.1750E+14 

.390 

28 

PW601 BOT 





1458 

1423 

1 372 



0020 







.044 

.1750E+14 

400 

29 

PW603 TOP 





1265 

1218 

3 310 



0005 







.036 

2470E4-14 

,780 

30 

PW603 MID 





1408 

1368 

3 310 



0010 







.042 

.2400E-fl4 

760 

31 

PW603 BOT 





1388* 

1338 

3 310 



0006 







.045 

.2470E-fl4 

.780 

32 

ORNL612 TOP 



5 

1240 

1143 

5 205 

.00290 

0005 







.043 

.4290E4-14 

2.540 

33 

ORNL612 MID 





1350 

1248 

5 250 



0037 









.4290E4-14 

2.540 

34 

ORNL612 BOT 





1285 

1188 

5 205 



0009 









.4290E^!4 

2.540 

35 

ORNL642 





1475 

1453 

9 966 



0034 









.lOSOE+14 

1.160 

36 

ORNL649 





1460 

1433 

10 512 



.0004 









. 1240E<-14 

1.430 

37 

ORNL652 MID 





1543 

1478 

3 180 



.0009 









.1960E>14 

.700 

38 

ORNL656 TOP 





1538 





.0003 









.1600Bn4 

.600 

39 

ORNL656 MID 





1538 





.0011 









.leOOE-vH 

.600 

40 

0RNL656 BOT 





1538 





.0013 









I600E-*-14 

.600 

41 

ORNL658 MID 





1480 

1458 

9 779 



.0014 









.9820E4-13 

1.080 

42 

ORNL658 BOT 





1460 

1438 

9 779 



.0009 









.9820E4-13 

1.080 

43 

ORNL660 MID 





1558 

1478 

2 635 



0005 









.2320E-fl4 

.070 

44 

ORNL662 TOP 





1363 

1338 

10 270 



.0001 









llSOEfU 

1.160 

45 

ORNL662 MID 





1498 

1473 

10 270 



0008 









1130Efl4 

1.160 

46 

ORNL662 BOT 





1488 

1463 

10 270 



0004 









.llSOEfM 

1.160 

47 

ORNL665 TOP 





1228 

1133 

9 583 



0136 









.4220Efl4 

4.580 

48 

ORNL665 MID 





1332 

1238 

9 583 



.0596 









.4220Efl4 

4.580 

49 

ORNL665 BOT 





1275 

1178 

9 583 



0178 









.4220Efl4 

4.580 

50 

ORNL667 MID 





1473 . 

1448 

6 209 



0002 









1090E-^14 

780 

51 

ORNL667 BOT 





1478 

1453 

6 209 



.0003 









.1090E+14 

.780 

52 

ORNL669 TOP 





1160 

1108 

10 352 



.0005 









.2330E*I4 

2.120 

53 

ORKL669 MID 





1300 

1253 

10 352 



.0013 









.2330E>14 

2.720 

54 

ORNL669 BOT 





1255 

1203 

10 352 



.0012 









.2SiOE*14 

2.720 

55 

ORNL002 3FMT 





1497 

1476 

10 173 



.0006 









.9820En3 

1 120 

56 

ORNL003 TOP 





1443 

1373 

11 985 



.0600 









.3090EfI4 

4.170 

57 

ORNL003 2FMT 





1503 

1433 




.1030 











58 

ORNL003 3FMT 





1528 

1458 




.1150 











59 

ORNL003 BOT 





1503 

1433 




.1260 
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ORNL UNI TOP 



. 7 

1393 . 

1313 

8 250 

..00293 

.0014 

.9271 

.7747 

.2667 

.060 

.18J0E-fl4 

1.780 

61 

ORNL UN3 TOP 



7 

1598 

1573 

5 807 

.00293 

.0010 

.9271 

.7747 

.2667 

.060 

.2350E4-14 

1.450 

62 

ORNL UN3 BOT 



3 

1920 

1673 

5 807 

.00293 

.0710 

.9271 

.7747 

0 

.055 

.2SS0E4-14 

1.750 

63 

NASA 504B 





1255 

1255 

8 070 

.00300 

.0005 

9110 

.7950 

.2550 

.060 

.5050E«13 

.470 

64 

NASA 503B 





1270 

1255 

8 070 

.00300 

0010 

.4580 

.3980 

.1420 

051 

.9560E^13 

.890 

65 

NASA 503C 

L 


L 


1271 

1255 

8 070 

.00300 

0010 

.4580 

.3980 

.1420 

.051 

.9840E^13 

.920 


®FueI code 2 refers to UN with axvcen content over 1000 ppm. 
*^Clad code 4: Nb-lZr 

CUd code 5: PWC-II (Kb-l wt. ^tZr-0 1 %C) 

Clad code 7: T-Ill (Ta-8 wt. ^^W-2 vx. 

Clad code 3 W-26Re (or W-25Re) 




TABLE X. - MATERIALS PROPERTIES 


Material 

Property 

Source 

Elastic modulus 

UN 

E = 1 38.7x10® - 100. 5x10® ^Vl - 6. 16xlO‘®(T - 77)] 

Refs. 20 and 21 


\ 


Nb-lZr 

E = 11. 8x10® - 3667 T 

Ref. 22 

PWC-11 

Assumed same as Nb-lZr 


T-111 

E = 26. IxlO® - 1700 T 

Ref. 23 

W-26Re 

E = 62.8X10® - 4667 T 

Private communication from I. Fiero, 



NASA Lewis Research Center 

Poisson s ratio 

1 



UN 

y = 0.284 - 0.382 -i- 

Ref. 20 


Vf 


T-111 

y = 0.25 

Unpublished TRW data obtained from 



I. Fiero, NASA Lewis Research Center 

Others 

y = 0.30 



-1 

Creep rate constants €^(hr ) = a^e 


UN 

e = 3.88X10-“ e-®10®®/T ^4.661 

e e 

Ref. 24 

Nb-lZr 

e = 3.88X10-“ e-'^4 604/T ^5.988 

e e 

Ref. 25 

PWC-ll 

c = 5.717x10-1® e-42WT ^5.988 

Private communication from R. L. Smith, 



NASA Lewis Research Center 

T-111 

e = 2.214X10-® e-'l'^ (7®-®®^ 
e e 

Ref. 25 

W-26Re 

e = 5.95xl0-®e-4®®''®/Ta®-®'^'' 

Private communication from W. F. Mattson, 



NASA Lewis Research Center 
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APPENDIX F 


SYMBOLS 


A creep rate constant 

88 ^ 

Aj^ activity of Kr in sweep gas, disintegrations /cm** -sec 

3 

si activity/10 

a^ constants (i = 1, 2, . . . ,7) 

B total birth (atoms formed by fission) , atoms 
b bxurnup, percent viranium atoms fissioned 

O 

c fission gas concentration, atoms/cm 

D total diffusion coefficient, cm /sec 

Djq fission-enhanced diffusion constant, cm^/fission 

2 

D atomic diffusion constant, cm /sec 

o ’ ' 

2 

D surface diffusion constant, cm /sec 

s ’ ' 

d combination of terms; 0. 359 fYV D7r^/4y| 

dp pore diameter, cm 

E elastic modulus 

F constant, fraction of fission gas released from surface 
^ axial force 

Q 

f fission-rate density, fissions/cm -sec 

f' f/10^^ 

g sum defined by eq. (Cl) 

h constant 

K gas constant, 1, 381x10“^^ N-cm/molecule-K 

Kj thermal conductivity, W/cm-K 

stress proportionality constants (i = 1,2, 3) 
constant 

L recoil length in fuel, cm 

Ng number of grains 
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n creep rate equation stress exponent 
P pressure 

p fraction of recoils stopped by gas after ejection 

energy of activation for atomic diffusion, J/mole 

Q energy of activation for surface diffusion, J/mole 
s 

R gas constant, 8. 314 J/mole-K 

Rg release rate from grain, atoms/sec 

R^p , release rate from interconnected porosity, atoms/sec 

Rpjj^ release rate from a fuel pin, atoms/sec 

Rg release rate from the fuel surface, atoms/sec 

R,p total gas release, atoms 

rj fuel inside radius, cm 

rg fuel outside radius or clad inside radius, cm 

rg clad outside radius, cm 

2 

S fuel surface area, cm 

T temperature , K 

O 

^ temperature/10 

t time , sec 

t^ decay time in sweep gas, sec 

u function of t 

3 

Vj volume of fuel, cm 

3 

V volume of intergranular gas, cm 

o 

3 

V initial volume of intergranular gas, cm 

V fuel bulk strain rate 

difference between calculated and measured values for i*’^ data set 

V pore velocity, cm/sec 

Jr 

V fission yield, atoms/fission 


y 

yi 


coordinate distance, cm 
thickness, cm 
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Y 2 half thickness, cm 

0 !j constants 

/3 proportionality constant, y = j3V /(V + V.) 

O & 

y fraction of intergranular voliime connected to fuel surface 
6 constant 

correction factors in Taylor expansion 
e creep rate 

X decay constant, sec“^ 

V Poisson's ratio 

7T dimensionless parameter 
<y stress 

T test time 

Subscripts : 
c clad 

e equivalent 

f fuel 

i internal 

o external 

r coordinate direction 

z coordinate direction 


9 


coordinate direction 
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